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ABSTRACT 

Aims. PKS 0548-322 (z = 0.069) is a "high-frequency-peaked" BL Lac object and a candidate very high energy (VHE, E > 100 GeV) y-ray 
emitter, due to its high X-ray and radio flux. Observations at the VHE band provide insights into the origin of very energetic particles present in 
this source and the radiation processes at work. 

Methods. We report observations made between October 2004 and January 2008 with the H.E.S.S. array, a four imaging atmospheric-Cherenkov 
telescopes. Contemporaneous UV and X-ray observations with the Swift satellite in November 2006 are also reported. 

Results. PKS 0548-322 is detected for the first time in the VHE band with H.E.S.S. We measure an excess of 216 y-rays corresponding to a 
significance of 5.6 standard deviations. The photon spectrum of the source is described by a power-law, with a photon index of F = 2.86±0.34stat ± 
O.lOsys. The integral flux above 200 GeV is ~ 1.3% of the flux of the Crab Nebula, and is consistent with being constant in time. Contemporaneous 
SwiftfKRT observations reveal an X-ray flux between 2 and lOkeV of F2_iokev = (2.3 + 0.2) x 10"" ergcm"^ s"', an intermediate intensity state 
with respect to previous observations. The spectral energy distribution can be reproduced using a simple one-zone synchrotron self Compton 
model, with parameters similar those observed for other sources of this type. 

Key words, gamma rays: observations - galaxies: active - galaxies: BL Lacertae objects: individual: PKS 0548-322 



1. Introduction 

BL Lacertae objects (BL Lacs) are an extreme class of Active Galactic 
Nuclei (AGN). The spectrum of these peculiar objects, which is flat 
and associated with compact sources in the radio, extends up to the 
y-ray band. None or weak emission lines are detected, and the radio 
and optical emission is highly polarized. BL Lacs are characterized by 
a rapid variability in all energy ranges, and display jets with appar- 
ent superluminal motions. The extreme properties of BL Lacs are ex- 
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plained by relativistic beaming, i.e. of a relati vistic bulk motion of the 
emitting region toward s the observer (see e.g. lBlandford & Reeslll97^ : 
lUrrv & Padovanill 1 9951) . 

The observed spectral energy distribution (SED) of BL Lacs usu- 
ally shows (in a vFy representation) two spectral components. The first 
peak is located in the radio to X-ray range, whereas the second is 
at higher energies, sometimes in the VHE range. The SED is com- 
monly explained by two different types of models. In leptonic mod- 
els, the lower energy peak is produced by synchrotron emission of 
relativistic leptons in a jet that points towards the observer. The sec- 
ond peak originates in the "inverse-Compton" scattering of leptons off 
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seed piiotons l lGhisellini & Maraschilll989l) . Depending on the origin 
of the seed photons, the leptonic models are divided in two classes. 
In the "synchrotron-self Compton" (SSC) models, the seed photon s 
come from the synchrotron photon field itself l lMarscher & Gear! 19851) . 
In external Compton (EC) scenarios, the seed photons are provided 
by various s ources, including the accretio n disk and broad emission 
line regions jPermer & Schlickeised Il993h . In hadronic models, the 
VHE emissi on is produced via the inte ractions of relativi stic protons 
with matter jPohl & Schlickeis ei||20Q(]| ), ambient photons ('Mannheim! 
119931) . magnetic fields ( I AharonianI 120001 : Piicke & Prot heroe 2001). 
BL Lac objects are divided into classes defined by the energy of the 
synchrotron peak: "low-energy-peaked" BL Lacs (LBLs) have their 
peak in the IR/optical wavelength whereas "high-energy-peaked" BL 
Lacs (HBLs) peak in the U V/X-ray band i fGiom mi & Padovani 1994: 

iPadovani & Giommill 19951) . 

PKS 0548-322 is a nearby (z = 0.069. iFosburv & Disnevlll976l) 
and bright BL Lac object hosted by a giant elliptical galaxy of abso- 
lute v isual magnitude My = -23.4 jFalomo et aLlll995l : IWurtz et alj 
[1993) which is the dominant member of a rich cluster of galaxies. The 
synchrotron emission of this obj ect peaks in the X-r ay band and there- 
fore it is classified as an HBL dPadovani & Giorrim i 1995). It is the 
third BL Lac detected by the HEAO X-ray satellite HMushotzkv et alj 
Il978h . Since then, PKS 0548-322 has been extensively studied by 
different X-ray experiments and satellites, showing a complex spec- 
tral behaviour. The X -ray spectrum deviates from a simple power- 
law dUrry et aLlll986h , and broken powe r-law models indica te a syn- 
chrotron peak energy in the 1-5 keV range. iTueller et aLl l l2008l) reported 
the detection of PKS 0548-322 in the hard X-ray band by Swif tfBAJ. 
Absorption features reported bv lSambruna & Mushotzkvl ( Il998h . were 
interpreted as the presence of circumnuclear ionized gas , but this has 
not b een confirmed by other spectroscopic observations l lBlustin et alj 
120041) . Thermal emiss ion from the host g alaxy on the kpc scale was 
detected by Chandra jDonato et al.ll2003l) . M ore details about the X- 
ray history of this object can be found in iPerri et aTl ( l2007h . The 
blazar was not detected in the MeV-GeV range by the EGRET de- 
tector jHartman et al.lll999l). PKS 0548-322 is sugg ested as a candi- 
date VHE emitter bv lCostamante & Ghisellinil ( l2002h . Previous obser- 
vations failed to detect y-ray emission from this object. Atmospheric 
Cherenkov telescope (ACT) experiments set upper limits on the VHE 
y-ray flux of this sou rce (CANGAROO: Fe>l5t,\ < 4.3x lO^'-cm-^s"', 
iRoberts et al.ll 999|: Durham Mark V I telescope: F£>3oogcV < 2.4 x 
10" cm -^s ' ■ IChadwick et al.l [2000l) . The discovery of VHE y-rays 
from PKS 0548-322 with the H.E.S.S. Cherenkov telescopes is pre- 
sented in this paper (see sections ICTI and \22\ . Contemporaneous obser- 
vations were carried out in X-ray, UV, and optical (see section lSTl 13.21 
and \33\ with the Swift satellite. In section|4]a single zone SSC model 
is applied to the data and discuss observational prospects with Fermi. 



2. H.E.S.S. Observations and Results 

2.1. H.E.S.S. observations 

The H.E.S.S. array of atmospheric-Cherenkov telescopes 
jAharonian et al.l l2006ij) observed PKS 0548-322 between October 
2004 and January 2008 with the full four-telescope array, for a total 
observation time of more than 63 hours. After the application of 
quality-selection criteria and a dead-time correction, a total of 34.9 
hours of good-quality data remain. The mean zenith angle of the 
observations is 10 degrees, which corresponds to a post-analysis energy 
threshold of 250 GeV. Th e results hereafter are based on t he H.E.S.S. 
3D-model analysis (see jLemoine^ Goumard et aP l2006h ). where a 
3D model is used to reconstruct the detected atmospheric shower 
induced by y-rays. For each detected shower, the direction, energy and 
3D-width are reconstructed. 

y-ray-like events are selected using cuts on image size, 3D-width, 
and telescope multiplicity. Only events that trigger at least three tele- 
scopes are kept in order to improve the y-hadron separation. For the 
analysis, on-source data are taken from a circular region of radius 9 
around the source position, and the background is subtracted using 
the event background rate estimated with the "ring background" model 
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Fig. 1. Distribution of squared angular distance ((f) for counts 
from for on-source and normalized off-source events. 



dAharonian et alj|2006ah . for which an annulus at the same distance to 
the center of the camera as the target position is used (excluding the 
region close to the source). The cuts applied for the analysis are sum- 
marised in Table [T] and are taken from lAharonian et" 'H] ( 120063) . 



Table 1. Summary of applied cuts. 



cut 


value 


image centroid distance 


2.5 deg 


image amplitude 


>60 photoelectrons 


# of telescopes 


>3 


core impact position 


<500m 


reduced 3D width a> 


<0.002 




<0.01 deg2 



2.2. H.E.S.S. Results 



6^) for counts from the 
1260 on-sources events 



The distribution of squared angular distance 
source is given in Figure [T] A total of A'on^ 
and A'ofF = 3105 off-source events are measured. The on-off normal- 
isation factor is a = non/f^off = 0.336, where SI is the solid an- 
gle of the respective on- and off-source regions. The observed ex- 
cess is Ny = A'on - tt^off = 216 y-rays, corresponding to a signifi- 
cance of 5.6 standard deviations according to Equation 6 from lLi & Mai 
| |1983|) . A two-dimensional Gaussian fit of the excess yields a position 
anm = 5''50'"38.4^±3.2\ feooo = -32° 16' 12.9"+40.2". The measured 
position is compatible with the nominal position of PKS 0548-322 
(aj20oo = 5''50"40.5%(5j2ooo = -32° 16' 16") at the 1 cr level. The 
differential photon spectru m shown in Figu re |2] is obtained using a 
forward-folding technique dPiron et al. II2000I) . Note that the likelihood 
maximisation provides the best set of parameters corresponding to a 
power law, with the coiTesponding error matrix. This is represented in 
Figure [2] as a 68% confidence-level band (grey zone). Points are de- 
rived from the residuals in different energy bins for illustration pur- 
poses only and plotted with 68% confidence-level error bars in binned 
flux . For a simple power-law hypothesis of the form E^^, a likelihood 
maximization yields a spectral index of 2.86 ± 0.34siai + O.lOsys and 
an integral flux: I(> 250 GeV )= (2.7 ± 0.6,m ± 0.5sys) x 10"'- cm"^ s"' 
iX-ld.o.f. = 24. 2/14). This corresponds to -1.3% of the H.E.S.S. Crab 
Nebula flux (see lAharonian et alj2006al) above the same threshold, well 
below the upp er limits reported by CANGAROO and the Durham Mark 
VI telescope dRoberts et al.ll 19991 : IChadwick et al.ll2000l) . No evidence 
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Fig. 2. VHE spectrum of PKS 0548-322. The shaded region rep- 
resents the Icr confidence level of the fitted spectrum, using a 
power-law hypothesis. 

for flux variability is seen in the data: fitting a constant flux to the data, 
rebinned on a yearly timescale, gives ax' of 2.1/3 d.o.f. Furthermore, 
there was no 1 hour detection > 5cr, which would be detected by 
H.E.S.S. if the source flux would have increased by a factor 4. Thus, 
a factor 4 of variability of the source flux during more than 1 hour is 
ruled out, and a flux variability factor of 1 + 3/ 's/t is excluded during 
a time T expressed in hours. 

3. Multi-Wavelength Analysis and Results 

3.1. Swift/XRT data reduction 

Obse rvations of PKS 0548-322 were carried out by Swift jGehrels et alj 
l2004h contemporaneously with H.E.S.S. for an effective 4.3 ks expo- 
sure on November 28, 2006 (obsid 0030836001). The XRT observatio n 
was performed in the Photon Count readout mode jHill et al.ll2004h . 
The XRT event files were calibrated and cleaned with standard filter- 
ing criteria with the xrtpipeline task using the latest calibration files 
available in the Swift CALDB distributed by HEASARC (SWIFT-XRT- 
CALDB-11 released in May 2008). Events in the ener gy range 0.3- 
lOkeV with grades 0-12 was used in the analysis (see lBurrows et all 
l2005l for a definition of XRT event grades). The source count rate was 
high enough to cause some photon pile-up in the inner 5 pixels ( ~ 9" 
) radius circle within the peak of the telescope Point Spread Function 
(PSF). Pile-up effect is avoided by selecting events within an annular 
region with an inner radius of 5 pixels and an outer radius of ~ 20 
pixels. The background is extracted from a nearby source-free circu- 
lar region of 20 pixel radius. Ancillary response files for the spec- 
tral analysis are generated with the xrtmkarf task, applying corrections 
for PSF losses (psfftag=yes). The PC mode grade 0-12 response ma- 
trix swxpc0tol2s0.20010101v01 1 .nnf is used in the fits. The spectrum 
is binned to ensure a minimum of 50 counts per bin and the 0.3-10 
keV recommended energy band is used. The flux is consistent with 
being constant during the observations. The spectral analysis is per- 
formed for the entire duration of the observation. A single power-law 
(F = 1.87 ± 0.06), with a Galactic absorption o f Nf = 2. 52 x 10^" cm^^ 
consistent with the Galactic value (Murphy e t alJll996h . gives a rea- 
sonable fit ix^ lA.o.f. = 40/25). The integrated energy flux between 
2 and lOkeV is Fa-iokeV = (2.3 ± 0.2) x 10"" erg cm"- s"', an in- 
termediate intensity state with respect to previous X-ray observations 
jBlustinet all 120041) . A broken power-law model with Galactic ab- 
sorption yields a relatively smooth spectral break (Fi = 1.7 ± 0.1, 
F2 = 2.0 ± 0.2, Zibrcak = 1-7;^Q 5 keV). However, it is not a significantly 
better fit (^^/d.o.f. = 35/25). In both cases, i.e. with a power-law or a 
broken power-law, the possible presence of an absorption feature near 
0.7 keV is observed. An F-test shows the addition of an absorption 
edged is significant at the 69% and 92% level for the power-law and 
broken-power-law cases, respectively. The presence of such an absorp- 
tion feature is discussed in a recent paper by B lustin et all j2004h . If real, 
it could be transient, as the analysis of XMM/RGS data by these authors 



reject the presence of such a feature in their spectra. In the present case, 
the significance is relatively weak. 

Although variability in optical, within a factor 2, has been reported 
for this source iPerri et al. 2007), no clear indication of X-ray variability 
has been found in Swift observations, confirming Perri et al.'s findings. 
This would indicate that the source has been observed in a relatively 
steady state, although the existence of more active periods can not of 
course be excluded 



3.2. Swift/UVOT data reduction 

The UltraViolet/Optical Telescope (UVOT) observations were made us- 
ing 5 filter settings, during the same time slots as the XRT data. Using 
the standard pipeline products (CALDB 20060917), the sequence of 
V, B, U, UVWl, UVM2 exposures have been checked for variability. 
No indication for variability is found and the individual frames are co- 
added. The reference stars are steady and UV OT photometry is w ithin 
0.03mag consistent with the literature values l lPoole et^l2008l) . The 
chosen aperture is 3" radius for all filters to avoid contamination by 
foreground stars and decrease the background noise. An aperture cor- 
rection!] estimated to be of the order of 10%, has been added to the data. 
A similar sized region is used to estimate the background. The result- 
ing count rate is corrected for Galactic absorption using a reddening of 
El% = 0.035 dSchlegel et al.lll998l) . consistent l lZombeck Il2007l) with 

the A'^''' from the power law fit to the Swi/Z/XRT data. The color in- 
dex is then convert ed into an extinction c oefficient using the conversion 
factor provided bv lGiommi et al ] <i2ooa> . The extinctions are shown in 
Tablell] 



Table 2. Extinction coefficients for PKS 0548-322 in the UVOT 
photometric bands. 



Filter 


A{A) 


A.i (mag) 


V 


5460 


0.112 


B 


4350 


0.147 


U 


3450 


0.182 


UVWl 


2600 


0.2345 


UVM2 


2200 


0.3395 


UVW2 


1930 


0.2905 



3.3. ATOM data reduction 

Optical observations were taken using the ATOM tele scope at the 
H.E.S.S. site from November 2006 until January 2009 jHauser et al.l 
.20041). Absolute flux values have been calculated using di fferential pho- 
tometry a gainst two stars calibrated bv lSmith et alj ( fl99ll) (B band) and 
IXie et al.1 Il996h (V,R,I band). A 4" radius aperture was used for all 
filter bands. 

A total of 312 measurements in 4 filter bands were taken in 
November 2006, May 2007, September 2007 to March 2008, and 
October 2008 to January 2009. Within errors, the photometry is compat- 
ible with a constant flux value of mg = 17.0, my = 16.2, niR = 15.6, and 
m, = 14.9 (4.4/6.7/8.1/10.4 x lO-'^ ergcm-^ s"') for the B, V, R and 
I band. It is compatible with the assumption that the optical emission 
comes dominantly from the host galaxy. B and V band are systemati- 
cally 0.2 mag brighter with respect to UVOT value s, but we consider it 
consistent within the instruments systematic errors l lPoole et aU2008l) . 



' http://heasarc.nasa.gov/docs/swift/analysis/threads/uvot_thread_ape 
rture.html 
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Table 3. Parameters used to fit the data (see text for further details). 
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Fig. 3. Spectral energy distribution of PKS 0548-322. Open circles: Radio, IR, and optical archival data are from the NED database. 
Fill red squares: ATOM data. Fill green circles: Swift observations contemporaneous with H.E.S.S. observations (corrected for 
Galactic absorption). Open diamonds: upper limits from SwiftfBALT, EGRET, Mark VI, and CANGAROO. The shaded region 
represents the Icr confidence bounds of the spectrum observed by H.E.S.S. above 250 GeV. The dotted lines at VHE show the 
spectrum of the source de-absorbed of the EBL attenuation. The left-hand side inlay details a port ion of the observ ed UV to X-ray 
spectrum. The dotted lines are the low level and high level flux observed with Swift and reported in lPerri etal ]l2007l The ri ght-hand 
side inlay details the H.E.S.S. spectrum. The upper limits deduced from first year Fermi observations ( Abdo et ani2009l) are also 
reported. 



4. Spectral Energy Distribution and discussion 

Figure [3] shows the SED of PKS 0548-322. Radio, IR, and optical 
archival data are taken from the NEEQ database. X-ray data from 
5w!/f/XRT are shown corrected for Galactic absorption, using the 
energy dependence for the photoelectric cross-section provided by 
XSPEC. The Swift/WOT data are corrected for Galactic absorption 
using the extinction coefficients i n Table [2] The Swift/ BAT flux level in 
the 14-195 keV range reported bv lTueller et all ( l2008h is interpreted as 
an upper limit in the SED, and this for several reasons. First, this flux is 
higher by a factor of 2 compared to the 2-10 keV XRT flux, and is not 
consistent with the spectral index after the energy break found in the 
SwiftpiRT analysis. Furthermore, these hard X-ray observations were 



^ The NASA/IPAC Extragalactic Database (NED) is operated by the 
Jet Propulsion Laboratory, California Institute of Technology, under 
contract with the National Aeronautics and Space Administration. 



integrated over 9 consecutive months that are not contemporaneous with 
the H.E.S. S. observations, and variability cannot be ruled out. Secondly, 
iTueller eF al. (2008) use a spectral index of 2.15 for all their sources. 
Using the BeppoSAX fit jPerri et alj|2007t) . the index would be 2.3-2.8 
in the energy range 14—195 keV. Hence, the harder assumption made by 
ITueller et alj ( 120081) probably ove restimates the actual fl ux in that band. 
The upper lirnits fr om EGRET dHartman et aDll999h. CANGAROO 
dRoberts et alj|199^ and Mark VI Durham telescope l lChadwick et all 
|2000|) are also displayed. The Icr confidence level of the VHE spectra 
measured by H.E.S.S. is shown, as well as the one year Fermi upper 
limits (Abdo et al. 200i), estimated at the decorrelation energy (where 
the upper limit doesn't depend on the assumed spectral index). 

A simple one-zone, homogeneous, time independent, synchrotron 
self-Compton (SSC) leptonic model is applied to interpret the contem- 
poraneous observations by H.E.S.S., Swift, and ATOM. The radio data 
is assumed to originate from an extended region farther along the jet and 
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the bump in IR-optical is attributed to emission from the host galaxy. 
Hence, these data are not fitted by the model. 

The SSC model describes the system as a spherical emitting re- 
gion of radius R, filled with a tangled magnetic field B. and prop- 
agating with a b ulk Dop pler factor d (e.g . [Band & Grindlav 1985; 
iKatarzvns ki et al. 2001: Gieb els et al]|2007l) . The energy distribution 
function (EDF) of the radiating leptons is described by a broken power- 
law, with indices «i and «2. between Lorentz factors y,,,;,, and y„iax, with 
a break energy ybreak and a normalisation factor K. The absorption by 
the extra-galacti c background light (EBL) is taken into account, adopt- 
ing the model of lFranceschini etaP d2008h . The luminosity distance is 
di = 9.6 X 10-^ cm, assuming the following cosmological parameters; 
Ho = 70kms-'Mpc-',A = 0.7,n = 0.3,qo = -0.5. 

The model parameters are underconstrained, and many different 
solutions provide a satisfactory fit of both X-ray and VHE spec- 
trum. We fit the data assuming that the flux did not exceed the (non- 
contemporaneous) Fenniupper limit at the time the H.E.S.S. observa- 
tions occured. The parameters used in the model are summarised in 
Table[3] Data can be fitted with a reasonable bulk Doppler factor <5 = 10. 
A lower value would result in a significant internal 7-7 absorption. 
The EDF is a broken power-law with indices «i = 2.2 and «2 = 3, 
constrained by the shape of the synchrotron spectrum. These values 
are compatible with the index of the broken power-law that fits the 
X-ray data from SwiftPiKT. We emphasize that the value of the break 
T2 - Fi observed by Swift is compatible (within errors) with 0.5, and the 
dominance of the synchrotron component compared to the Compton 
one is consistent with a relativistic plasma losing its energy primarily 
via synchrotron cooling. The EDF is defined between 7n,i„ = 1 and 
7max = 5 X 10*. The value of the maximal Lorentz factor y^^,, is not 
well constrained due to the lack of contemporaneous observation in the 
hard X-ray range. The break energy in the EDF is associated with the 
synchrotron break that is detected around 4 x 10" Hz. In the model, the 
energy break depends on the magnetic field of the emitting region with 
the following relation; ybieak (6(5)"''", and is thus well constrained for 
a given Doppler factor. By choosing a value for the magnetic field B, 
the two last parameters of the model (A^ and R) are constrained by the 
flux in X-ray and at VHE. 

With the adopted parameters, the emitting region is out of equipar- 
tition with a ratio of the particle energy density to the magnetic en- 
ergy density Ue/us = 30. However this high value is in agreement with 
that usually derived in VHE blazars using one-zone models. One can 
also compute from these parameters a synchrotron cooling timescale 
fcooi ~ 280 ks at the energy break, smaller than the source crossing time 
fcross = R/c = 640 ks. Both values must be divided hy 6 = 10 to give 
an apparent dynamical timescale t^.^ = R/Sc ~ 64 ks, and an apparent 
cooling timescale ~ 28 ks , both being less than one day. The absence 
of a clearly detected short term variability does not allow to put a strong 
constraint on these parameters, alt hough historical ob servations do ex- 
hibit some variability at year scale l lBlustin et alj2004l ). which is clearly 
compatible with these values. The radius of the emitting source derived 
from the model is about a few hundreds times the gravitational radius 
of the blazar's cen tral black hole, assuming a mass of Mbh ~ 10**^ Mo 
jBarth et ani2003l) . A larger Doppler factor could give an equally good 
fit with a smaller radius, a smaller break Lorentz factor, a larger B field, 
but a higher ratio of particle to magnetic energy density, so more out 
of equipartition. However these solutions would be favored if shorter 
variability timescales were detected. In any case, the low Compton to 
synchrotron flux ratio indicates that the source is (on average ) far from 
the Compton catastrophe, although more intense, short-lived flares can 
not be excluded. 

5. Conclusion 

Observations performed by H.E.S.S. from 2004 up to 2008 have es- 
tablished PKS 0548-322 as a VHE 7-ray source amongst the closest 
blazars. The VHE flux is consistent with being constant within the 
H.E.S.S. observation period. The contemporaneous X-ray data can be 
fit by a broken power-law model with Galactic absorption. 
For the first time for this source, a SED comprising contemporaneous 
optical, UV, X-ray and VHE measurements is made. A one-zone SSC 



model, taking into account absorption by the EBL, provides a satisfac- 
tory description of these data. Observations of this object with Fermi 
should help to place stronger constraints on the model parameters by 
building a more complete SED of the object. It is a further confirmation 
that the proportion of TeV sources among close blazars is fairly high. 
The Doppler factor used in the modeling by a simple one-blob model 
is compatible with the value commonly adopted for relativistic jets, al- 
though, as is often the case for TeV sources, no superluminal motion 
has been reported for this object. The addition of new objects to the set 
of "high-frequency-peaked" BL Lacs observed at VHE, combined with 
Fermi data, should soon enable population studies and insights into the 
underlying physical processes. 
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